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Abstract  We  have  addressed  the  derivation  and  the  analysis  of  an  adaptive 
decision  scheme  to  detect  possible  extended  targets  modeled  as  Gaussian  vectors 
known  to  belong  to  a  given  subspace;  noise  returns  from  the  cells  under  test 
are  modeled  as  independent  and  identically-distributed  Gaussian  vectors  with 
one  and  the  same  covariance  matrix;  a  set  of  secondary  data,  free  of  signal 
components  is  also  available;  secondary  data  are  Gaussian-distributed  and  share 
the  same  covariance  matrix  of  noise  in  the  cells  under  test  but  for  a  possible 
different  power  level. 

The  proposed  detector  relies  on  a  two-step  design  procedure:  first  we  derive  the 
GLRT  assuming  that  the  noise  covariance  matrix  is  known  up  to  a  scale  factor; 
then,  we  come  up  with  a  fully  adaptive  detector  by  replacing  the  structure  of  the 
covariance  matrix  of  the  noise  with  the  sample  covariance  matrix  based  upon  the 
secondary  data.  The  first  step  requires  the  maximum  likelihood  (ML)  estimate  of 
the  covariance  matrix  of  the  useful  signal  (under  the  signal-plus-noise  hypothesis) 
which,  in  turn,  has  a  known  structure.  That  ML  estimate  has  been  firstly  proposed 
by  Bresler  in  [3];  a  different  derivation  is  also  proposed  herein. 

The  performance  assessment  is  conducted  resorting  to  the  method  proposed  in 
[4-5]  to  model  extended  targets:  therein  an  exponential  model  for  fully-polarized 
returns  has  been  used  assuming  that  each  scattering  center  can  be  characterized 
by  its  (relative)  range,  amplitude,  and  polarization  elipse. 
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ABSTRACT 

This  paper  addresses  adaptive  detection  of  extended  tar¬ 
gets  whose  dominant  scatterers  are  modeled  as  a  linear 
combination  of  known  modes  with  jointly  Gaussian  co¬ 
efficients.  At  the  design  stage  we  consider  the  following 
instances  of  a  priori  knowledge  about  target  and  noise: 
the  noise  is  white  with  unknown  power  level  while  the  co- 
variance  matrix  of  the  coefficients  is  either  known  up  to 
its  eigenvalues  or  totally  unknown  (for  the  first  case  we 
also  assume  orthogonal  modes);  or  both  the  noise  and 
the  coefficients  covariance  matrices  are  unknown.  We 
derive  GLRT-based  detectors  for  the  first  two  cases  and 
propose  a  fully-adaptive  detector  for  the  third  one  by 
resorting  to  secondary  data  possessing  the  same  covari¬ 
ance  matrix  of  the  data  under  test  up  to  a  multiplicative 
factor.  Finally,  we  assess  the  performance  of  the  last  two 
detectors  by  resorting  to  Monte  Carlo  simulation. 

1.  INTRODUCTION 

A  High-Resolution  Radar  (HRR)  can  resolve  a  target  into 
a  number  of  scattering  centers  depending  on  the  range  ex¬ 
tent  of  the  target,  the  range  resolution  capabilities  of  the 
radar,  and  its  operating  frequency.  Measurements  indicate 
that  radar  properties  of  several  targets  can  be  modeled  in 
terms  of  a  set  of  scattering  centers  each  parameterized  by 
its  range,  amplitude  and,  possibly,  polarization  elipse  [1], 

Properly  designed  HRRs  allow  significant  enhancement 
of  the  detection  performance  as  shown  in  [2,  and  references 
therein]  where  adaptive  radar  detection  of  extended  deter¬ 
ministic  targets,  embedded  in  possibly  non-Gaussian  distur¬ 
bance,  has  been  addressed.  Therein,  possible  returns  from 
target’s  scattering  centers  are  modeled  as  signals  known  up 
to  multiplicative  factors,  possibly  different  from  one  scat- 
terering  center  to  another,  namely  supposed  to  belong  to 
a  one-dimensional  subspace  of  the  observables.  From  this 
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viewpoint,  a  more  general  framework  is  considered  in  [3] 
where  adaptive  detection  of  point-like  deterministic  targets, 
assumed  to  belong  to  a  known  subspace  of  the  observables, 
has  been  addressed.  That  model  has  also  been  generalized 
to  deal  with  point-like  stochastic  targets  in  [4].  Herein,  fol¬ 
lowing  the  lead  of  [2,  4],  we  address  adaptive  detection  of 
extended  targets  modeled  as  Gaussian  signals  known  to  be¬ 
long  to  a  known  subspace  of  the  observables. 

To  be  more  definite,  assume  that  an  array  of  L  anten¬ 
nas  senses  I\  range  cells  and  that  each  antenna  collects  T 
samples  from  each  of  those  cells;  thus,  denoting  by  X;,.  the 
A'-dimensional  vector,  with  N  =  LT,  containing  returns 
from  the  k- th  cell,  k  =  1, . . . ,  K,  the  problem  at  hand  is 
that  of  detecting  the  possible  presence  of  the  target’s  scat¬ 
tering  centers  s&  G  C'V,  k  =  1, ...  ,K,  modeled  as  a  lin¬ 
ear  combination  of  r  linearly  independent  modes,  namely 

JV”  x  f 

as  s j.  =  H0fc,  where  the  columns  of  the  matrix  H  G  <L 
are  the  vectors  of  the  basis  while  the  entries  of  the  vector 
Ok  G  C  are  the  coordinates  of  the  possible  target’s  scat- 
terer  in  the  /.'-th  cell.  Moreover,  we  assume  that  the  0ks 
are  (complex)  Gaussian  vectors  with  zero-mean  and  positive 
semi-definite  covariance  matrix  Rgg  =  E[0k0 jj,  where  E 
denotes  statistical  expectation  and  t  Hermitian  transpose, 
i.e., 

Ok  ~  CJ\fr(0,  Reg), 

and  that  the  0k\  are  independent.  The  noise  vectors  ri/,  G 
C\  k  =  1  ,...,K,  are  modeled  as  additive  and  indepen¬ 
dent  complex  Gaussian  vectors  with  zero-mean  and  covari- 

•  •  i  T-»  /r^NxN  „ 

ance  matrix  given  by  7Rnn  G  (L.  ,  7  >  0,  i.e.  ~ 

CA/7v(0, 7R,m). 

In  this  paper  we  design  detectors  based  upon  the  Gener¬ 
alized  Likelihood  Ratio  Test  (GLRT)  assuming  the  follow¬ 
ing  instances  for  target  and  noise: 

•  Case  1:  Rm,  =  IjV,  where  Ijy  denotes  the  N  x  N 
identity  matrix,  7  is  unknown,  and  the  covariance  ma¬ 
trix  Rgg  of  the  0fcS  is  known  up  to  its  eigenvalues, 
ei, . . . ,  er  say.  We  also  assume  that  the  columns  of 


H  are  each  other  orthogonal. 

•  Case  2:  Rnn  =  I  .y  while  7  and  the  covariance  matrix 
Ree  are  unknown. 


matrix,  E  €  Crx'  is  a  diagonal  matrix  with  non-negative 
diagonal  entries,  ei, . . . ,  er,  and  1  denotes  conjugate  trans¬ 
pose. 

It  follows  that 


•  Case  3:  both  7R nn  and  R00  are  unknown,  but  M  > 
N  zero-mean  Gaussian  vectors  with  covariance  ma¬ 
trix  R„n  are  available.  These  vectors  will  be  referred 
to  in  the  following  as  secondary  data. 


where 


Rss  =  HEH*. 


H  =  HU  = 


hi 


eC 


Nxr 


We  also  assess  the  performance  of  detectors  designed  to 
tackle  cases  2  and  3  by  resorting  to  Monte  Carlo  simula¬ 
tion.  The  derivation  of  the  GLRT  for  case  2  (unknown  R^e) 
is  conducted  taking  advantage  of  the  results  reported  in  [5]. 

More  precisely,  the  next  section  addresses  the  design  of 
the  GLRT-based  detectors  while  Section  3  contains  a  pre¬ 
liminary  performance  assessment  conducted  by  resorting  to 
Monte  Carlo  simulation. 


is  a  slice  of  a  unitary  matrix  (remember  that  HH  =  Ir). 
Moreover,  if  we  denote  by  H-1  G  ^<nx(n  r)  ^  matrjx 
whose  columns  are  orthonormal  vectors  spanning  the  or¬ 
thogonal  complement  of  the  space  described  by  the  columns 
of  H,  the  following  identity  holds  true 


+  7Ijv 


HH1 


E  +  7I?- 

7liV-r 


Ht 

H^t 


2.  DETECTOR  DESIGNS 


It  follows  that  the  pdf  of  X,  under  Hi,  is  given  by 


We  cannot  implement  the  likelihood  ratio  since  total  or  par¬ 
tial  ignorance  of  R.on,  7,  and,  possibly,  R„„  is  assumed. 
Leaving  aside  case  3,  we  can  resort,  instead,  to  the  GLRT 
which  is  tantamount  to  replacing,  under  each  hypothesis, 
the  unknown  parameters  with  the  maximum  likelihood  esti¬ 
mates. 

To  this  end,  note  that  the  detection  problems  to  be  solved 
can  be  re-cast  in  terms  of  the  following  binary  hypothesis 
test: 

(H0:  Xk~CATN  (0,7ljv),.  k  =  1, . . . ,  K, 

[H1  :  xfc  ~  CM N  (0,RSS  +7Iat)  ,  k=l,...,K, 

where 

Rss  =  E[sfesj.]  =  HReeH^. 


/(X|ffll7)E) 


-1 _ 1  ir _ I _ 

nNK  ry(N  r)K  11  _|_  ei)K 


(2) 


x  e 


-iy,f  x;pixfc-yK  yr 

7  Z— lk  =  l  kT,  =  1 


where  P~  and  P~  denote  the  projectors  onto  the  subspace 

spanned  by  h,  and  the  orthogonal  complement  of  the  sub¬ 
space  spanned  by  the  columns  of  H,  respectively. 

For  further  processing,  it  is  also  convenient  to  denote  by 
i  =  1, , . .  ,r,  the  function 

*. 

1  _  k  h» 

5i(X  |  7,e»)  =  - - -ff  e  ? 

(7  +  ei) 


2.1.  Case  1:  7  and  the  as  unknown 

This  subsection  is  aimed  at  designing  a  GLRT  to  deal  with 
the  case  that  7  and  the  eigenvalues  of  R.oo  are  unknown 
(while  the  eigenvectors  of  R.00  are  known). 

For  further  developments  it  is  necessary  to  determine  the 
probability  density  function  (pdf)  of  the  data  matrix,  i.e., 

X  =  [xi  •  •  'Xjf], 


under  each  hypothesis.  The  pdf  of  X  under  Hq  is  given  by 


/(X|ff„,7) 


1 

^NK 


NK 


(i) 


The  pdf  of  the  data  under  the  II\  hypothesis,  instead, 
can  be  determined  following  the  lead  of  [4].  To  this  end, 
denote  by  UEU^  an  eigendecomposition  of  the  positive 
semidefinite  matrix  R00;  hence,  U  €  Crxr  is  a  unitary 


and  re-cast  the  pdf  of  the  data  matrix,  under  Hi,  as 


f(X\HU7,V)  = 


X 


1  1 

^NK  j(N-r)K  e 


1  X^K 

1 


Pixfc 

H 


n^(xl  7,  ei). 

i= 1 


The  derivation  of  the  GLRT  requires  substituting  the  un¬ 
known  parameters  by  the  corresponding  maximum  likeli¬ 
hood  estimates  under  each  hypothesis.  Maximizing  the  pdf 
of  the  data  under  H0  with  respect  to  7  is  straightforward;  it 
yields 


max  f(X  |  H0, 7) 


1  _NK  (  NK 

NK  \  XfXl 

\Z^k=  1 


NK 


Now  focus  on  the  more  difficult  task  to  maximize  the 
pdf  of  the  data  matrix  with  respect  to  the  e^s,  i  =  1, . . . ,  r, 


and  7,  under  the  Hi  hypothesis.  We  attack  this  problem  by 
maximizing  the  pdf  with  respect  to  the  e.,s,  given  7.  To  this 
end,  it  is  easy  to  check  that 


max  ft  (X  |  7,  e-i)  = 

ei>  0 


e  if  7  <  lir 


7 1 

1  _ 

zice  * 


if  7  >  7 i, 


where 


k= 1 


K 


Moreover,  note  that 


i  =  1 . r. 


maxft(X  |  7, e,),  *  =  l,...,r, 

ei>0 

is  constant  for  7  <  7^,  while  its  restriction  to  [7$,  +00 [  is  a 
monotonically  decreasing  function. 

Now,  denote  by  7(7  the  i-th  order  statistic  obtained  by 
ranking  71, . . . ,  7r  in  ascending  order,  i.e., 


7(i)  <■■■<!( 


Moreover,  let 


K 


7  = 


ExipgXfc- 


(N  —  r)K  ^ 

Then,  it  is  easy  to  check  that  the  function 


Proof.  Let 

MX  |  7)  = 


1  1  K(N-r)-, 

re  7 


7 tNK  7(JV-r)X' 


n 


—  K 


K 


n 


(7(0 ) 


1  K^(i) 

Ke  7 


j  =  1, . . .  ,r  +  1,  with 


TT  1  1 

llyfe  7  =1 

ie0  ' 


and 


n< 


—K 


(7  w) 


K 


=  1. 


It  is  easy  to  check  that 

maxei>0,...,eT.>o  /(X  |  fTi,7,  E) 

f  /ii(X|7),  if  0  <  7  <  7(i), 

MX  |  7),  if  7(i)  <  7  <  7(2), 


=  < 


(4) 


MX  I  7),  if  7(r—  1)  <  7  <  7(r), 
l  Mi(x  |  7),  if  7  >  7(r)- 


max  /(X|J?1>7,E) 

ei>0,...,er>0 

attains  its  maximum  (with  respect  to  7)  at  7  if  7  <  7m, 
while  the  maximizer  is  a  point  of  the  interval  [7(1) ,  7]  if  7  > 
7(7).  More  specifically,  it  is  possible  to  prove  the  following 
theorem. 

Theorem  1.  The  function 


max  /(X  |  #7,7,1}),  7  >  0, 

ei>0,...,er>0 

attains  its  maximum  with  respect  to  7  at 


7(1), 

if 

7(1)  <  7(i), 

7(0’ 

if 

7(<)  >  7(0,*  =  !,  •  •  •  J  ~  ^%)  <  7 (j). 

7(r+l) 

,  if 

7(0  >7(0,*  =  1,---,L 

where 

{N- 

r)7  +  J2i=i  7(0  .  J  r  +  1 

-r  +  (j-l)  ’  J  1, . . .  ,r  +  f. 

7(0  = 

N 

Moreover,  the  hj(X  |  7)  possess  the  following  properties: 

•  h:j  ( X  |  7)  >  /i)+7  (X  |  7),  7  >  0,  with  equality  if 

7  =  7 (j)J  =  1,  •  ••  ,r; 

•  hj(X.  |  7),  7  >  0,  attains  its  absolute  maximum  at 

7  (j)J  =  !,•••,  r  +  1. 

It  follows  that 

•  if  7(1)  <  7(i)  *en 

max  /(X  |  i?i,  7,  E),  7  >  0, 

ei>0,...,er>0 

attains  its  maximum  at  7(7)  =  7. 

•  If  7(7)  >  7(i)  and  7(2)  <  7(2),  it  is  also  true  that 

7(i)  <  7(2)  <  7(2); 

moreover,  since 

max  hi(X  |  7)  =  h±(X  |  7(1))  =  h2(X  |  7(1)) 

7<7(l) 

<  max  /i2(X  |  7) 

7(l)<7<7(2) 

=  h2(X  |  7(2)), 


we  conclude  that 


where 


max  /(X  |  . . .  ,er),  7  >  0, 

ei>0,...,er>0 

attains  its  maximum  at  7 (2).  Moreover,  for  j  =  3, . . . ,  r, 

if  7(i)  >  7(i)>  i  =  1,  •  •  •  ,3  -  1.  and  7(j)  <  j(j),  it  is 

also  true  that 


yu- 1)  <  %)  <  7(jj; 


moreover,  since 


max  /ii(X  |  7) 

7<7(l) 


=  /ii(X  I  7(1))  =  fe2(X  I  7(1)) 
<  max  /i2(X  |  7) 

7(l)<7<7(2) 

=  h2(X  |  7(2)) 


^  *4  vw 

x  /  N-r  +  (j- 1)  \ 

V^-r^+ECiSw/ 


(N-r+U-l))K 


and  7 (j)  is  given  by  eqn.  (3). 


2.1.1.  GLRT-based  detector 

It  is  then  easy  to  see  that  the  decision  statistic  of  the  GLRT 
is  given  by 

A[X] 


we  conclude  that 


< 

r  At  [x] , 

7(i) 

< 

7(i), 

A2[X], 

7(i) 

> 

7(1), 7(2)  < 

7(2), 

= 

hj{X 

1  7(j-i))  =< 

A  j  [X] , 

%) 

> 

7 (»),*  =  1,  • 

•  j  - 

VI 

< 

max  hj  (X  |  7) 

7(r-l)<7<7(i) 

= 

hj(X 

l7(i))- 

Ar  [X] , 

%) 

> 

7 (»),*  =  1,  • 

..,r- 

'T(r) 

7T 

VI 

,  1, 

7 (i) 

> 

7  (»),*  =  1,- 

. .  ,r, 

where 


max  /(X  |  i?r,7,E),  7  >  0, 

ei>0,...,er.>0 

attains  its  maximum  at  J(j). 

Similarly,  if  7V)  >  77  j,  i  =  1, ... ,  r,  we  conclude 
that 


Aj[X]  = 


(me  Efc=ixlxfe 

j- 1 


(AT-r)7+J] 


7  (i) 


Af-r+(j-l) 


(JV-r+0-t))if 


nr.,  (7.,)) 


max  /(X  |  ^,7,11),  7  >  0, 

ei>0,...,eT.>0 

attains  its  maximum  at  7(r+1). 

The  theorem  is  thus  proved.  ■ 

It  follows  that 


2.2.  Case  2:  7  and  Rw  unknown 

Again  the  pdf  of  the  data  matrix  X  under  //(J  is  given  by 
eqn.  (1).  It  follows  that  implementation  of  the  GLRT  for  the 
case  at  hand  requires  determining  the  pdf  of  X  under  H 1 
and  maximizing  it  with  respect  to  7  and  .  The  pdf  of  X 
under  II  \  is  given  by 


max  max  f(X  I  iTi  ,7,  E), 

7>0  ei>0,...,er>0  v  1  '  ’ 


is  obtained  by  evaluating  the  function  (4)  at  its  maximizer; 
we  get 

max7>o  max  /(X|ffi,7,E) 
el  >0,...,er>0 


'  /i(X|Jfi),  7(1)  <  7(1), _ 

/2(X|iCi),  7(1)  >  7(1), 7(2)  <  7(2), 


=  {  /j(X|/fi), 


7 «)  >  7(i) ,  *  =  1,  •  ■  ■ ,  J  —  1, 7(j)  <7 (j), 


/r(X|J7l),  7(i)  >  7(i),»  =  ■  ■  ,r  -  l,7(r)  <  7(r), 

l  /r+i(X|Ri),  7(q  >  7(i) ,  *  =  1,  •  ■  ■ ,  r, 


/(X|ff1)7,Ree) 


,,-NK 


|R  SS  +  7Jiv 


if 


e-A'tr[R(R„+7IN)-1] 


where  tr[-]  and  |  •  |  denote  the  trace  and  the  determinant  of  a 
square  matrix,  respectively,  while 

R  =  — XXf. 

K 

The  maximization  of  the  above  pdf  with  respect  to  7  >  0 
and  Ree  >  0  has  been  addressed  in  [5]  where  it  has  been 
shown  that 


max  max  /(X|fTi,  7, Roe)  =  maxmax g(X|fTi,  7,  Q) 

7>0  R0<9>O  7>0  Q>0 


where  g(X.\Hi,  7,  Q)  is  given  by 


It  follows  that  the  statistic  of  the  GLRT  is  given  by 


0(X|ff1)7,Q) 


g-ftrpp^R+fQ+i,.)-1?] 

[(7r7nQ  +  Ir|f 


with 

Q  =  L^ReeL, 

7 

P  =  L_1HtRHL_t, 

and  L  denoting  a  nonsingular  r  x  r  square  root  factor  of 
HtH,  namely  such  that  LL0  =  H0H.  Moreover,  maxi¬ 
mization  of  g(X.\Hi,  7,  Q)  has  been  performed  in  two  steps: 
first  it  is  shown  that 


maxQ>0 


5(X|ff1)7)Q) 

5r(X|7),  7  —  5 

5r-l(X|  7),  0r  <  7  < 


5l(X|7),  02  <  7  <  01, 

,  5o  (X I  7),  7  >  0i, 


where 


5j(X|7) 


e~jK 

(irj)NK 


0-f{tr[PiR]+E’'=J+1 


and  ,(f>r  are  the  eigenvalues  of  the  matrix  P  ordered 

in  non-ascending  order,  i.e., 


f  A  r  [X] ,  7r  :  0r; 

Ar— l[x],  7r  7  0n  7^— 1  —  0r— 15 


A[X]  =  { 


where 


Ai[X],  7^  >  0j,i  =  2, . . .  ,r,7r  <  0i, 

1,  7i  >  0i,i  =  1,.  ..,r, 


(5) 


AjfXl  = 


(  jvV  Sfe=i  x!xfc 


NK 


This  detector  will  be  referred  to  as  a  partially-adaptive  one. 


2.3.  Case  3:  7,  Rnn  and  R^e  unknown 

As  a  final  step  observe  that  the  above  decision  statistic  can 
be  modified  to  address  detection  in  presence  of  correlated 
noise.  To  this  end,  assume  that  M  >  N  additional  Gaussian 
vectors  17  ~  CM  n  (0,  Rrm),  k  =  K  +  1, . . . ,  K  +  M,  free 
of  signals  components,  but  possessing  the  same  statistical 
characterization  of  the  noise  in  the  cells  under  test,  but  for  a 
possible  different  power  level,  are  available.  These  vectors 
are  usually  referred  to  as  secondary  data.  Now  we  have  to 
discriminate  between  the  Hq  hypothesis  that  can  be  re-cast 
as 

(rk  ~  CM n  (0)  7R  nn)  1  k  =  1, . . . ,  K, 

\r fc  ~  CM n  (0,  Rrm)  j  k  =  K  +  1, . . . ,  K  +  M, 

and  the  Hi  hypothesis  which  is  given  by 


0r  <  •  •  •  <  01. 

Then,  it  is  shown  that  the  following  theorem  holds  true. 
Theorem  2.  The  function 

max 5(X|Ri,7,Q),  7  >  0, 

attains  its  maximum  with  respect  to  7  at 


Vk  ~  CM N  (0,  R$s  T  7Rnn)  ,  k  —  1,  .  .  .  ,  if, 

Yk  ~  CM N  (0i  Rrm)  ,  k  =  K  +  1,..  .  ,K  +  M, 

where  R,m  denotes  the  unknown  covariance  matrix  of  the 
secondary  data  and,  up  to  a  multiplicative  factor,  that  of  the 
data  under  test. 

Then,  we  can  come  up  with  a  fully-adaptive  detector 
by  resorting  to  detector  (5)  fed  by  the  whitened  data  x*,  = 
Rnn'  Yk,  where  Rnn  is  the  sample  covariance  matrix  based 
upon  secondary  vectors,  i.e.. 


7 ri  'Ir  —  0r? 

7r— 1,  7 r  7  0r,7r— 1  ^  0r— 1, 

7i>  7i  >  0i.*  =  2,...,r,7i  <  0i, 
,  7o,  7*  >  <t>i,i= 


where 


r 

i=j  + 1 


0i  +  tr[PjjR] 

N-j 


j  =  0,  •  •  - ,  r. 


j  k+m 

=  M[  rferfc' 

k=K+ 1 

This  detector  will  be  referred  to  as  a  fully-adaptive  detector. 

2.4.  Performance  Assessment 

In  the  following  we  assess  the  performance  of  the  partially- 
adaptive  and  the  fully-adaptive  detectors  by  resorting  to  stan¬ 
dard  Monte  Carlo  counting  techniques.  More  precisely,  in 
order  to  evaluate  the  threshold  necessary  to  ensure  a  pre¬ 
assigned  value  of  Probability  of  False  Alarm  (Pfa)  and  the 


Proof.  See  [5],  ■ 


Figure  1 :  Pd  vs  SNR  (dB)  of  the  partially-adaptive  detector 
for  N  =  16,  r  =  3,  Pfa  =  10-3,  7=1 ,  K  and  pg  as 
parameters:  curves  with  a  square  at  each  data  point  refer  to 
pg  =  0.5  while  those  with  a  star  correspond  to  pg  =  0.9; 
solid,  dashdot,  and  dotted  refer  to  K  =  24,  K  =  12,  and 
K  =  6,  respectively. 


Figure  2:  Pd  vs  SNR  (dB)  of  the  partially-adaptive  and 
the  fully-adaptive  detector,  N  =  16,  r  =  3,  pg  =  0.5, 
pn  =  0.6,  7=1,  Pfa  =  10-3,  K  as  parameter:  curves  with 
a  square  at  each  data  point  refer  to  the  partially-adaptive 
detector  while  those  with  a  star  correspond  to  the  fully- 
adaptive  one  with  M  =  32;  solid,  dashdot,  and  dotted  refer 
to  K  =  24,  K  =  12,  and  K  =  6,  respectively. 


Probability  of  Detection  (Pd)  we  resort  to  100 / Pfa  and  1000 
independent  trials,  respectively.  It  is  important  to  stress  that 
both  detectors  have  the  Constant  False  Alarm  Rate  (CFAR) 
property  with  respect  to  7. 

Figures  1  and  2  plot  Pd  vs  the  signal-to-noise  ratio  (SNR) 
defined  as 


tr(7R„„) 

with  Rfm  =  Ijv  for  curves  of  Figure  1 .  More  precisely,  the 
curves  of  Figure  1  refer  to  the  partially-adaptive  detector 
with  IV  =  16,  r  =  3,  7  =  1,  an  exponentially-correlated 
target  with  one-lag  correlation  coefficient  pg,  i.e., 

R -ee(i,j)  =  OgplT'71,  G  {1, 

K  and  pg  as  parameters.  As  to  H  its  columns  are  the  first 
r  vectors  of  the  canonical  basis  in  R^.  The  figure  high¬ 
lights  that  the  performance  improves  as  pg  increases:  as  a 
matter  of  fact,  since  the  noise  is  white  the  two  hypotheses 
become  more  and  more  distinguishable  as  the  signal  corre¬ 
lation  increases.  In  addition,  the  performance  of  the  detector 
improves  as  K  increases.  In  Figure  2,  instead,  we  plot  the 
performance  of  the  fully-adaptive  detector  implemented  by 
resorting  to  M  =  32  secondary  data  and  assuming  7=1 
and  an  exponentially-correlated  noise  with  one-lag  correla¬ 
tion  coefficient  pn,  i.e., 

=  Pn~J\  *>J  G  {1 

For  comparison  purposes  curves  of  the  partially-adaptive 
detector,  fed  by  whitened  data,  are  plotted  too.  The  fig¬ 
ure  highlights  the  loss  of  the  fully-adaptive  detector  with  re¬ 
spect  to  the  partially-adaptive  one  where  an  unknown  noise 


covariance  must  be  estimated  from  secondary  data.  A  thor¬ 
ough  performance  assessment  is  the  object  of  outgoing  re¬ 
search  activity. 

3.  REFERENCES 

[1]  W.  M.  SteedlyandR.  L.  Moses,  “High  Resolution  Ex¬ 
ponential  Modeling  of  Fully  Polarized  Radar  Returns,” 
IEEE  Trans,  on  Aerospace  and  Electronics  Systems, 
Vol.  27,  No.  3,  pp.  459-469,  May  1991. 

[2]  E.  Conte,  A.  De  Maio,  and  G.  Ricci,  “CFAR  de¬ 
tection  of  Distributed  Targets  in  non-Gaussian  distur¬ 
bance,”  IEEE  Trans,  on  Aerospace  and  Electronic  Sys¬ 
tems,  Vol.  38,  No.  2,  pp.  612-621,  April  2002. 

[3]  S.  Kraut,  L.  L.  Scharf,  and  L.  T.  McWhorther,  “Adap¬ 
tive  Subspace  Detectors,”  IEEE  Trans,  on  Signal  Pro¬ 
cessing,  Vol.  49,  No.  1,  pp.  1-16,  January  2001. 

[4]  L.  T.  McWhorter  and  L.  L.  Scharf,  “Matched  Subspace 
Detectors  for  Stochastic  Signals,”  11th  MIT  Lincoln 
Labs  Workshop  on  Adaptive  Sensor  and  Array  Process¬ 
ing,  Lexington,  MA  (USA),  March  11-13,  2003. 

[5]  Y.  Bresler,  “Maximum  Likelihood  Estimation  of  a  Lin¬ 
early  Structured  Covariance  with  Application  to  An¬ 
tenna  Array  Processing,  Fourth  Annual  ASSP  Workshop 
on  Spectrum  Estimation  and  Modeling,  Minneapolis, 
MN  (USA)  pp.  172-175,  August  3-5,  1988. 


Abstract 


Presentation 


Back  to  TOC 
Next  Abstract 


Adaptive  Radar  Detection  of  Extended  Gaussian  Targets 


by 


Giuseppe  Ricci  and  Louis  L.  Scharf 


Problem  Formulation 


This  paper  addresses  adaptive  detection  of  extended  Gaussian  vectors 
known  to  belong  to  a  known  subspace  of  the  observables. 

To  this  end,  we  assume  that: 

•  an  array  of  L  antennas  senses  K  range  cells  and  that 

•  each  antenna  collects  T  samples  from  each  of  those  cells. 

Denote  by  the  LT-dimensional  vector  containing  returns  from  the  k- th 
cell,  k  =  1  The  problem  herein  addressed  is  the  detection 

of  the  possible  presence  of  the  target’s  scattering  centers  s^  E  C  ,  k  = 
1 modeled  as  a  linear  combination  of  r  linearly  independent  modes, 
namely  as 

within  the  K  range  cells  under  test. 

We  assume  that: 

-  the  matrix  H  E  CNxr  is  known. 

-  The  coordinates  vectors  Ok  ~  £A/*r(0,  R^)  are  each  other  independent 
and  R00  is  a  positive  semidehnite  covariance  matrix. 

-  The  noise  vectors  n&  ~  CA/*at(0,  yRnn)  are  each  other  independent. 


Problem  Formulation  (cont.) 


We  design  detectors  based  upon  the  Generalized  Likelihood  Ratio  Test 
(GLRT)  assuming  the  following  instances  for  target  and  noise: 

•  Case  1:  we  assume  that 

-  the  scale  factor  7  is  unknown; 

-  H  is  a  slice  of  a  unitary  matrix  and  is  known; 

-  the  covariance  matrix  of  the  Ok s  is  known  up  to  its  eigenvalues, 
ei, . . . ,  er  say. 

•  Case  2:  we  assume  that 

-  7R nn  =  7l^V; 

-  the  scale  factor  7  is  unknown; 

-  the  covariance  matrix  R^^  is  unknown. 

•  Case  3:  we  assume  that 

both  7Rnn  and  R00  are  unknown,  but  M  >  N  zero-mean  Gaussian 
vectors  with  covariance  matrix  Rnn  are  available.  These  vectors  will 
be  referred  to  in  the  following  as  secondary  data. 


Detector  Design:  7  and  the  e^s  unknown 


The  detection  problems  to  be  solved  is 

I  H0  :  Xfc  ~  CM  n  (0, 7ljv) ,  k  =  l,...,K, 

\  Hi  :  Xfc  ~  CM n  (0,  +  7ljv)  ,  k  =  1, . . . ,  I\, 

where  7  and  the  eigenvalues  of  R00  are  unknown. 


It  is  necessary  to  introduce  the  following  quantities: 

•  UEtA  is  an  eigendecomposition  of  R#,?; 

•  H  =  HU=  [hi-  hr]; 

•  7p)  is  the  i-th  order  statistic  obtained  by  ranking  the 

1  k 

7i  =  vE4p  h,x%  ^  = 

k= 1 

in  ascending  order  where  Pg  is  the  projector  onto  the  subspace  de¬ 
scribed  by  hp 

•  7  and  7(j)  are  given  by 


7  = 


1 

{N  -  r)K 


E 

k=  1 


and 


{N  ~  r)7  +  eLi  7(i) 
TV  —  r  +  (j  —  1) 


j  =  l,...,r  +  1, 


respectively.  As  to  Pj~  it  denotes  the  projector  onto  the  orthogonal 
complement  of  the  subspace  spannned  by  the  columns  of  H. 


Detector  Design:  7  and  the  e^s  unknown  (cont.) 


It  is  possible  to  prove  that  the  GLRT  is  given  by 


A[xi,...,xfc]  <  A 
Ho 


where 


A[xi,  •  •  •  ,xfc; 


with 


Ai[xb  . . .  ,Xjfc], 
A2[xi,  . . .  ,Xjfc], 

<  Aj  [xi , . . . ,  Xfc] , 

Ar[xi,  •  •  .  ,Xjfc], 

1, 


7(1)  <  7(1)7 

7(i)  >  7(i)  ?  7(2)  <  7(2)7 


7 (i)  >7p),^  =  l,---,i-l,7(i)  <7(j) 


7(*)  >  7(7  i  =  i,  •  •  ■ ,  r  -  1, 7(r)  <  7(r) 
7(i)  >  7(») 7  *  =  l,  -  -  - , 


Aj[xi, . . .  ,Xfe 


iVv  ^=1 4xfc 


NK 


/  (Af-r)7+ELi7(»)^ 
Y  iV-(r-j+l)  ^ 


(TV— (r— j+l))A' 


nr  • 


Detector  Design:  7  and  unknown 


The  detection  problems  to  be  solved  is 

Ho  :  Xfc  ~  CM n  (0, 7R) , 

H\  :  Xfc  ~  CM n  (0,  HR^R  +  7ljv)  , 

with  7  and  R00  unknown. 

It  is  necessary  to  introduce  the  following  quantities: 

•  R  is  given  by 

_  1  K 

R  =  77  E  **4; 

H  k= 1 

•  L  is  a  nonsingular  r  x  r  square  root  factor  of  RR  namely  such  that 

Lid  =  RH; 

•  the  r  x  r  matrix  P  is  dehned  as 

P  =  L^RRHL-t; 

•  (f)\, ...  Mr  denote  the  eigenvalues  of  P  ordered  in  non-ascending  order. 


k  =  1, . . . ,  K, 
&  =  !>■■■>  K, 


Detector  Design:  7  and  unknown  (cont.) 


It  is  possible  to  prove  that  the  GLRT  is  given  by 


A[xi  •  •  -Xjfc]  <  A 
Ho 


where 


'  Ar[xi  •  •  -Xfc], 
Ar-i[xi  •  •  -xfc], 


A[xi  •  •  •  xfc]  =  < 


Ai[xi  •  •  -xfc], 


7r  ^ 

7 r  7  7r— 1  —  0r— 1? 

7*  >  <f>i,i  =  2,  •  • ,  ,r,7i  <  0i, 
7*  >  (/>i,i  =  1,  •  •  •  ,r, 


where 


and 


r  1  Wv  vt.  ^ 


Aj  [xi  •  •  -xfc]  = 


yNK  ^ k= 


Hk=i  x[.xfc  J 


iN-j)KnL 


I,  = 


Vr 

^i=j+l 


+  trfPjjR 


)K 

l 


N-j 


j  =  1, 


r. 


This  detector  will  be  referred  to  as  a  partially-adaptive  one. 


Detector  Design:  7R nn  and  R^  unknown 


Assume  that  M  >  N  additional  Gaussian  vectors  CM  N  (0,  Rnn) , 

k  =  K  +  1, . . . ,  K  +  M,  are  available. 

The  detection  problem  to  be  solved  is 

rfc  ~  CM N  (0, 7 Rnn)  ,  k=  1,  •  •  • ,  K, 

*7  ~  CM N  (0,  R nn)  j  k  =  K  1,  .  .  .  ,  K  M, 

Yk  ~  CM  N  (0,  Rss  +  7Rnn)  ,  k=l,...:K: 

17  ~  CM N  (0,  R nn)  5  k  =  K  +  1, . . . ,  K  +  M, 

with  7Rnn  and  R00  unknown. 

We  can  come  up  with  a  fully-adaptive  detector  by  resorting  to  the  partially- 
adaptive  detector  (1)  fed  by  the  whitened  data  =  HM2rk,  where 

„  X  K+M 
Rnn  =  VT  X!  rfcrfc- 
k=K+ 1 


This  detector  will  be  referred  to  as  a  fully-adaptive  detector. 


Performance  Assessment 


Performance  assessment  is  conducted  via  Monte  Carlo  Simulation:  we  re¬ 
sort  to  105  and  103  independent  trials  in  order  to  set  the  threshold  necessary 
to  ensure  Pfa  =  10”3  and  Pd,  respectively. 

We  assume  the  following  simulation  parameters 

•  the  columns  of  H  are  the  first  r  vectors  of  the  canonical  basis  in 

•  the  target  is  exponentially-correlated  with  one-lag  correlation  coeffi¬ 
cient  pe,  i.e., 

I Mbi)  =  (Jept~3\  hi  e  {l,...,r}; 

•  we  choose  N  =  16,  r  =  3,  7  =  1; 

•  curves  of  Figure  1  assume  R nn  =  In  while  those  of  Figure  2  refer  to 
exponentially-correlated  noise  with  one-lag  correlation  coefficient  pnj 
i.e., 

R,™  '•/•  I' . Vj: 

•  the  fully- adaptive  detector  of  Figure  2  assumes  M  =  32  secondary 
data. 


Performance  Assessment  (cont.) 


Pd  vs  SNR  (dB)  of  the  partially-adaptive  detector  for  N  =  16,  r  =  3, 
Pfa  =  10-3,  7  =  1 ,  K  and  pe  as  parameters:  curves  with  a  square  at 
each  data  point  refer  to  pe  =  0.5  while  those  with  a  star  correspond  to 
pe  =  0.9;  solid,  dashdot,  and  dotted  refer  to  K  =  24,  K  =  12,  and  K  =  6, 
respectively. 


Performance  Assessment  (cont.) 


Pd  vs  SNR  (dB)  of  the  partially-adaptive  and  the  fully- adaptive  detector, 
N  =  16,  r  =  3,  po  =  0.5,  pn  =  0.6,  7=1,  Pfa  =  10-3,  K  as  parameter: 
curves  with  a  square  at  each  data  point  refer  to  the  partially-adaptive 
detector  while  those  with  a  star  correspond  to  the  fully-adaptive  one  with 
M  =  32;  solid,  dashdot,  and  dotted  refer  to  K  =  24,  K  =  12,  and  K  =  6, 
respectively. 
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Conclusions 


•  Note  that  both  the  partially-  and  the  fully-adaptive  detectors  have  the 
Constant  False  Alarm  Rate  (CFAR)  property  with  respect  to  7. 

•  A  thorough  performance  assessment  is  the  object  of  outgoing  research 
activity. 


